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The nature of the ordering of the one-dimensional Heisenberg spin-glass model with a long- 
range power-law interaction is studied by extensive Monte Carlo simulations, with particular 
attention to the issue of the spin-chirality decoupling/coupling. Large system sizes up to 
L = 4096 are studied. With varying the exponent cr describing the power-law interaction, 
we observe three distinct types of ordering regimes. For smaller cr, the spin and the chirality 
order at a common finite temperature with a common correlation-length exponent, exhibiting 
the standard spin-chirality coupling behavior. For intermediate cr, the chirality orders at a 
temperature higher than the spin, exhibiting the spin-chirality decoupling behavior. For 
larger cr, both the spin and the chirality order at zero temperature. We construct a phase 
diagram in the cr versus the temperature plane, and discuss implications of the results. 
Critical properties associated with both the chiral-glass and the spin-glass transitions are 
also determined. 
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1. Introduction 

In spite of quite extensive studies for years, the true nature of the ordering of spin-glass 
(SG) magnets still remains elusive and controversial. Since the magnetic interaction in most 
of real SG materials is known to be nearly isotropic, they should be described as a first 
approximation by the isotropic Heisenberg model. Recently, consensus appears among various 
numerical works that the isotropic Heisenberg SG in three dimensions (3D) exhibits a finite- 
temperature transition, while the nature of the transition still remains controversial.^"^) 

It has been suggested in Ref.^^ that the model might exhibit an intriguing "spin-chirality 
decoupling" phenomenon, i.e., the chirality exhibits the glass order at a temperature higher 
than the standard SG order, Tqg > Tsg-^^^^ Chirality is a multispin variable representing the 
handedness of the noncollinear or noncoplanar spin structures induced by frustration. Based 
on such a spin-chirality decoupling picture of the 3D Heisenberg SG, a chirality scenario of 
experimental SG transition has been advanced.^' ^' ^'^^ By contrast, Refs.^"^-* claim that the 
3D Heisenberg SG does not exhibit such a spin-chirahty decoupling, only a single transition 
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which is simultaneously SG and chiral-glass (CG). 

To get deeper insight into the behavior in physical dimension d = 3, it is often useful to 
study the phenomena by extending the dimensionality to general d dimensions. In the limit of 
low d, the short-range (SR) Heisenberg SG exhibits only a T = transition in d = 1. In d = 2, 
recent calculations suggest that the vector SG model, either the three-component Heisenberg 
SG^^) or the two-component XY SG/^^ exhibits a T = transition accompanied by the 
spin-chirality decoupling, i.e., the CG correlation-length exponent vcG is greater than the 
SG correlation- length exponent usg- The spin-chirality decoupling associated with a finite- 
temperature transition could occur, if any, in d > 3. In the opposite limit of high d, the SR 
Heisenberg SG model in infinite dimensions d — > oo reduces to the mean-field (MF) model, 
i.e., the Sherrington-Kirkpatrick (SK) model. The Heisenberg SK model is known to exhibit 
only a single finite-tcmpcrature SG transition, with no spin-chirality decoupling. In high but 
finite d, Monte Carlo (MC) result of Ref.^^^ suggested that the spin-chirality decoupling did 
not occur in d = 5, but might occur in d = 4. Reflecting an intrinsic difficulty in thermalizing 
large systems in high dimensions, however, the true situation still remains largely unclear. 

In the present paper, we attack the issue of the spin-chirality coupling/decoupling in the 
Heisenberg SG from a different perspective. Namely, we study a different type of Heisenberg 
SG model, i.e., the one-dimensional (ID) Heisenberg SG with a long-range (LR) power-law in- 
teraction proportional to (r is the spin distance). A preliminary account of the simulation 
was presented in Ref."'^*^-' 

In the limit of sufficiently large a — )• oo, the model reduces to the standard d = 1 model 
with a SR interaction. In the opposite limit of cr — >^ 0, the model reduces to an infinite-range 
model, i.e., the SK model corresponding to d = oo. (Note that for cr < 1/2, one needs to 
rescale the interaction strength by an appropriate power of L to make the energy extensive.) 
Hence, varying a of the ID LR model might be analogous to varying d in the SR model. 
Indeed, this correspondence was supported by recent studies by Katzgraber and Young^^' 
and by Leuzzi et al^^'^^^ for the Ising SG. Indeed, these authors have suggested more detailed 
correspondence between d of the SR model and a of the ID LR model, e.g., (i) the upper 
critical dimension d = 6 corresponds to o" = 2/3, (ii) the lower critical dimension, which lies 
between d = 2 and 3, corresponds to o" = 1, and (iii) d = 3 corresponds to cr ~ 0.9. 

Advantages of studying such ID models might be threefold. First, systems of large linear 
size L, never available in high dimensions, can be studied (up to L = 4096 in the present 
calculation). Second, one can continuously change and even fine-tune the parameter a playing 
the role of effective "dimensionality", while it is impossible to continuously change the real 
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dimensionality d in the SR model. Hence, by studying the properties of the ID model with 
varying a, one might get an overall picture concerning how the "coupling vs. decoupling" 
behavior depends on the effective dimensionality. Third, certain analytical results based on 
the renormalization-group (RG) calculations are available in ID, which might be utilized in 
interpreting the numerical data. 

Indeed, RG calculations, though did not take account of the possibility of the spin-chirality 
decoupling, suggested that the model exhibited a rich ordering behavior with varying a.^^'^^^ 
For (7 < 2/3, the Gaussian fixed point is stable and the model exhibited a finite-temperature 
SG transition of the MF type. For 2/3<cr<l, a non-trivial LR fixed point becomes stable 
leading to a finite-temperature SG transition characterized by the non-MF exponents. In 
particular, the critical-point-decay exponent is determined solely by the power describing the 
spin-spin interaction, and is given by rjsG = 3 — 2(T.^°'^^^ For a >1, the SG transition occurs 
only at zero-temperature with ijsg = 1- 

Meanwhile, it remains to be seen how the spin-chirality decoupling arises in this ID model 
with varying a. Since the MF Heisenberg SK model does not show the spin-chirality decou- 
pling, the spin-chirality decoupling associated with a finite-temperature transition should be 
realized, if any, only in the intermediate range of a, near or below a = 1. Thus, we study here 
both the spin and the chiral orderings of the model by large-scale MC simulations, varying a in 
the range 0.7 < a < 1.1, which spans the non-MF regime. By studying such a onc-dimcnsional 
Heisenberg SG model with LR power-law interactions, we are able to study the spin-chirality 
decoupling/coupling phenomena from a wider perspective. Our numerical results indicate that 
the model exhibits the spin-chirality decoupling in the range 0.8 a "S 1.1, while the usual 
spin-chirality coupling behavior occurs for a < 0.8. 

The paper is organized as follows. In §2, we define our model and explain some of the details 
of our numerical method employed. Various physical quantities calculated in our simulations 
are introduced in §3. Our MC results are presented in §4. Quantities like the CG and SG 
correlation-length ratios, the CG and SG susceptibilities, the CG and SG Binder ratios, the 
CG and SG overlap distribution functions, etc, are calculated for various values of the range 
parameter a. Then, a phase diagram of the ID LR model is constructed in the a versus 
the temperature plane. Critical properties associated with the CG and SG transitions are 
analyzed in §5 my means of a finite-size scaling analysis. Finally, §6 is devoted to summary 
and discussion. 
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2. The model and the method 

We study the ID classical Heisenberg model with the random LR power-law interaction 
Jij, whose Hamiltonian is given by 

H = ~ ^ JijSi ■ Sj , (1) 

<ij> 

where Si = {Sf,S^,Sf) is a three-component unit vector at the z-th site, and the < ij > 
sum is taken over all spin pairs on the chain once. The coupling Jij decays with a geometric 
distance rij as a power-law, 

•^ij = C = ^ — -p^, (2) 

'li Z^<ij>'ij 

where is an independent random Gaussian variable with zero mean and standard deviation 
unity. Periodic boundary condition is applied by placing L spins on a ring. Then, the geometric 
distance between the spins at i and j is given by rij = (L/tt) sin(7r \i — j\ / L). 

We perform extensive MC simulations for various values of a = 0.7, 0.8, 0.85, 0.9, 0.95, 
1.0 and 1.1. Preliminary result for a = 1.1 was reported in Ref.^^^ The lattice sizes studied 
are L= 128, 256, 512, 1024, 2048 and 4096. In our simulation, we use a single-spin-flip heat- 
bath and an over-relaxation method combined with temperature-exchange technique.^^^ We 
perform over-relaxation sweeps 5 times per every heat-bath sweep, which constitutes our unit 
MC step. 

Equilibration is checked by monitoring: i) All the "replicas" travel back and forth many 
times (typically more than 10 times) along the temperature axis during the temperature- 
exchange process between maximum and minimum temperature points, whereas the relaxation 
due to single-spin flip is fast enough (both chiral and spin autocorrelation times about 20 MC 
steps or less) at the maximum temperature: (ii) All the measured quantities converge to stable 
values. Error bars of physical quantities are estimated by the sample-to-sample statistical 
fluctuation over bond realizations. 

In Tables I and II, we show the details of our simulation conditions, including the system 
size L, the number of independent samples (bond realizations) Ng, the number of tempera- 
ture points used in the temperature-exchange process Nt, and the minimum and maximum 
temperatures and T^ax- In Table I, we show these conditions for a = 0.7, 0.8, 0.85, 0.9, 
and 0.95 which correspond to the Tsg > regime, and those for a = 1.0 and 1.1 in Table II 
which correspond to the Tsg = regime. 
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a 


L 


Ns 


Nt 


-'-max 


-'-mm 


0.7 


128 


896 


16 


0.275 


0.120 




256 


896 


16 


0.275 


0.120 




512 


896 


16 


0.275 


0.120 




1024 


512 


16 


0.275 


0.120 




2048 


256 


16 


0.275 


0.120 


0.8 


128 


896 


16 


0.190 


0.090 




256 


896 


16 


0.190 


0.090 




512 


896 


16 


0.190 


0.090 




1024 


512 


16 


0.190 


0.090 




2048 


256 


16 


0.190 


0.090 




4096 


256 


8 


0.190 


0.134 


0.85 


128 


896 


20 


0.170 


0.065 




256 


896 


20 


0.170 


0.065 




512 


896 


20 


0.170 


0.065 




1024 


512 


20 


0.170 


0.065 




2048 


512 


10 


0.170 


0.108 




4096 


256 


8 


0.170 


0.119 


0.9 


128 


896 


20 


0.150 


0.055 




256 


896 


20 


0.150 


0.055 




512 


896 


20 


0.150 


0.055 




1024 


896 


20 


0.150 


0.055 




2048 


896 


13 


0.135 


0.072 




4096 


256 


8 


0.135 


0.093 


0.95 


128 


896 


16 


0.113 


0.050 




256 


896 


16 


0.113 


0.050 




512 


896 


16 


0.113 


0.050 




1024 


896 


16 


0.113 


0.050 




2048 


512 


16 


0.113 


0.050 




4096 


256 


12 


0.113 


0.062 



Table I. Various parameters of our Monte Carlo simulations of the ID LR model with a = 0.7, 0.8, 
0.85, 0.9 and 0.95 which correspond to the Tsg > regime. L is the linear dimension of the 
system, Ng is the number of samples, Tmax and Tmin are the highest and the lowest temperatures 
used in the temperature-exchange run, and Nt is the total number of temperature points. 
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a 


L 


Ns 


Nt 


T 

^max 


-'-mm 


1.0 


128 


896 


16 


0.097 


0.0337 




256 


896 


16 


0.097 


0.0337 




512 


896 


16 


0.097 


0.0337 




1024 


896 


16 


0.097 


0.0337 




2048 


256 


32 


0.101 


0.0337 


1.1 


128 


512 


32 


1.0 


0.014 




256 


1024 


32 


1.0 


0.014 




512 


1024 


16 


0.079 


0.014 




1024 


448 


64 


0.22 


0.019 




2048 


256 


56 


0.22 


0.019 



Table II. Various parameters of our Monte Carlo simulations of the ID LR model with a = 1.0 
and 1.1 which correspond to the Tsg = regime. L is the linear dimension of the system, Ng is 
the number of samples, T^ax and T^in are the highest and the lowest temperatures used in the 
temperature-exchange run, and Nt is the total number of temperature points. 



3. Physical quantities 

In this section, we define various physical quantities calculated in the following section. 
The local chirality at the i-ih. site Xi is defined for three neighboring Heisenberg spins by 
the scalar 

Xi = Si+i ■ {Si X Si-i) . (3) 

First, we define an 'overlap' for the chirality. We prepare at each temperature two inde- 
pendent systems 1 and 2 described by the same Hamiltonian (1) with the same interaction set. 
We simulate these two replicas 1 and 2 in parallel with using different spin initial conditions 
and different sequences of random numbers. 

The /c-dependent chiral overlap, Qxik), is defined as an overlap variable between the two 
replicas 1 and 2 as a scalar 

1 ^ 

«xW = ^Exf^Xre*-, (4) 

i=l 

where the upper suffixes (1) and (2) denote the two replicas of the system, and is the 
distance along the chain (ring) in units of lattice spacing. 

The fc-dependent spin overlap, Qapik), is defined by a tensor variable between the a and 
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P components of the Heisenberg spin, 

1 ^ 

9"Mfc) = ]vE^a^i?«'"'^ («,/3 = x,2/,z). (5) 

1=1 

In term of the A;-dependent overlap, the CG and SG order parameters are defined by the 
second moment of the overlap at a wavevector k = 0, 

m _ [(kx(o)P)] 

ICG - =4 ' 



a,l3=x,y,z 



(2) 

The CG order parameter q^Q has been normalized here by the mean-square amplitude of the 
local chirality, 

1 ^ 

x' = 3]vEKxf)], (8) 

i 

which remains nonzero only when the spin has a noncoplanar structure locally. The local 
chirality amplitude depends weakly on the temperature and the lattice size, in contrast to the 
Heisenberg spin variable whose amplitude is fixed to be unity by definition. 

One often uses the Binder ratios to estimate the critical temperature. The CG and the 
SG Binder ratios are defined by 



^.c = ,(ll-9,^^). (10) 



lfii_9iM2)!)l 

2 V [{ism? 

These quantities are normalized so that, in the thermodynamic limit, they vanish in the 
high-temperature phase and gives unity in the non-degenerate ordered state. In the present 
Gaussian coupling model, the ground state is expected to be non-degenerate so that both gcG 
and gsG should be unity at T = 0. 

Finite-size correlation length of the ID LR model is defined by 

" 2sin(A;m/2j 

where the q{k)'^ is defined via cq.(4) for the chirality and by eqs.(5) and (7) for the spin, with 
kjn = The reason for the power l/{2a— 1) appearing in eq.(ll) is that, at long wavelength, 
we expect a modified Ornstein-Zernike form for the LR modeP^^ 



q^^\k)<x{t + k'"'-Y\ (12) 



^ 2sin(W2) U('?(fcm)2)] 7 ' ^^^^ 
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where t is a measure of the deviation from the critical point. 

The correlation length divided by the system size ^/L, the correlation-length ratio, is a 
dimensionless quantity. Around the critical temperature Tc, this quantity is expected to obey 
the finite-size scaling form, 

^ = X{{T-Tc)L^/n{l + aL-^), (13) 
ij 

where v is the correlation-length exponent, X a scaling function, a a constant, and u the 
leading correction-to-scaling exponent. 

We also calculate the CG and SG susceptibilities defined by 

XCG = Lq^cG ' ^SG = Lqfl- (14) 

While the SG and CG susceptibilities are dimensionfull quantities, they can be made 
dimensionless by dividing them by L'^~^ where 77 is a critical-point-decay exponent. Generally, 
the exponent rj is not known in advance, but in the case of the present LR interaction, rjsa is 
determined by the power describing the spin-spin interaction and is given by rjsa = S—2a?^''^^^ 
Thus, one expects the finite-size scaling form to hold around the critical temperature, 

= Y{{T - rc)LV-)(l + a'L-), (15) 

where F is a scaling function and a' is a constant. As L — >■ 00, the ^/L and the x/(L^~^) 
curves of different L plotted versus the temperature should asymptotically cross at T = T^. 
Unfortunately, exact expression is not known for the corresponding chiral-glass exponent r]cG- 
The chiral-overlap distribution P{qx) is defined by 

^'(i) = [(%;-9x(o)))]- (16) 

The spin-overlap distribution is defined originally in the tensor space with 3x3 = 9 compo- 
nents. To make this quantity more easily visible, one may define the diagonal spin-overlap, 
which is a trace of the original tensor overlap as^' 

P{qdiag) = mqdiag- ^ ^mmW))]- (17) 

IJ,=x,y,z 

4. Monte Carlo results 

In this section, we present the results of our MC simulation on the ID LR Heisenberg SG 
for various values of the range parameter a. 

4.1 a = 0.90 

We begin with the case of cr = 0.90, at which a the spin-chirality decoupling is observed 
most clearly. In Fig.l, we show the correlation-length ratios for the chirality ^cg/L (a), and 
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Fig. 1. The correlation-length ratio for the chirality (a) and for the spin (b) plotted versus the 
temperature for a = 0.9. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition 
point. Note that the S^cg/L data are given on a semi-logarithmic plot. 



0.7 




T 



Fig. 2. The spin-glass susceptibility ratio xsg/L'^~^^'^ plotted versus the temperature for a = 0.9, 
with an "exact" exponent value tjsg = 3 — 2cr = 1.2. The crossing points are indicated by small 
blue arrows. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. 




Fig. 3. The Binder ratio for the chirality (a) and for the spin (b) plotted versus the temperature for 
a = 0.9. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. 
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for the spin ^sg/L (b). As can be seen from the figure, the spin ^sg/L curves cross at 
progressively lower temperature as L increases, whereas the chiral ^cg/I^ curves intersect at 
an almost L-independent temperature. 

We plot in Fig. 2 the temperature dependence of the SG susceptibility ratio XSG / 
where 2 — rjsc = 2a — 1 = 0.8 for the present value of cr = 0.9. As can be seen from the figure, 
the data of different L exhibit a crossing behavior as expected for the dimensionless quantity 
with a finite-T transition. While the crossing points of isc/I-' approach the bulk SG transition 
point from above, those of XSG / approach the SG transition point from below. For the 




0.006 



0.008 



Fig. 4. The (inverse) size dependence of the crossing temperatures of ^cg/^, Csg/-^^7 and xsg/-Z^^ 

for a = 0.9. Lines represent power-law fits of the form (18). The CG and SG transition tempera- 
tures are extrapolated to be Tea = 0.105 ± 0.003 and Tsg ^ 0.086 ± 0.003. 




Fig. 5. The chiral-glass susceptibility ratio Xcg/ L^^'^^^ plotted versus the temperature for a = 0.9, 
where the value of rjcG is set to 1.9 as determined by the finite-size-scaling analysis of §5. The 
red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. The inset exhibits the 
(inverse) size dependence of the crossing temperatures in which the line represents a power-law fit 
of the form (18). The CG transition temperature is extrapolated to be Tcg = 0.103 ± 0.003 
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CG susceptibility ratio, on the other hand, this type of analysis has only restricted utility 
because of the lack of our knowledge of the chirality-chirality interaction and the associated 
r/cG-value. 

The Binder ratios are shown in Fig. 3 for the chirality (a), and for the spin (b). The 
chiral Binder ratio gcG exhibits a shallow negative dip, which tends to be even shallower 
with increasing L, in contrast to the behavior observed for qqq of the 3D SR Heisenberg SG 
where qcg of different sizes exhibited a crossing behavior on the negative side of gcG-^^ In 
contrast to this, gcG of the present ID model hardly exhibits a clear crossing behavior for 
smaller lattices. Interestingly, however, large lattices of L = 2048 and 4096 eventually exhibit 
two crossings, one at T ~ 0.11 on the negative side of gcG and the other at T ~ 0.09 on the 
positive side of gsG'- See the inset of Fig.3(a). 

In order to estimate the bulk CG and SG transition temperatures quantitatively, we 
need to extrapolate the crossing temperatures of either the correlation-length ratio or the SG 
susceptibility ratio to L = oo. In Fig.4, we plot Tcross{L) of the ^cg/L curves between of 
the sizes L and 2L as a function 1/L, together with the corresponding ones of the £^sg/L 
and the xsg/ L"^'^^^ curves. The error bar of each Tcross{L) is estimated on the basis of the 
bootstrap method over available samples, combined with a polynomial fit (of the fourth or 
the fifth order) of the temperature dependence of the physical quantities. 

An extrapolation to L = cx) is made on the basis of the relation, 

Tcross{L) — Tcrossioo) ^ cL ^ , (18) 

where c is a constant. For both cases of the SG and the CG, the exponent 9 is equal to 




Fig. 6. Overlap distribution function of the chirality (a) and of the spin (b) for a — 0.9 at a temper- 
ature T — 0.055 below Tsg- The sample average of i = 2048 is taken here for a subset of total 
samples (256 samples) 
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^ + w for the crossing temperature Tcross{L) of either the correlation- length ratio or the 
glass-susceptibility ratio. 

In the case of the SG, we perform a combined fit of both ^sg/L and of Xsg/ ^ 
where a common Tcrossipo) = Tga and Osg are assumed. We then get TgQ = 0.086ib 0.003 and 
0SG = 0.44ib0.07: See Fig.4. The errors shown here and below are obtained via the standard x^- 
analysis as explained in some detail in Ref.^^ The x^-value per degree of freedom of our fit here 
is x^/DOF 2± 1.24, with the associated fitting probability Q 2± 0.29, which is quite reasonable. 
The smallness of the error bar associated with Tsg comes from the fact that a combined fit of 
the two independent of Tcross{L), i.e., those of £,sg/L and of xsg/L"^'^^^ , is used here, each 
approaching Tsg either from above or from below. For the CG, we have Tcross{L) of ^cg/L 
only, since r]cG is not known in advance. The CG transition temperature is then estimated 
via a power-law fit of Tcross 

(L) of Ccg/L to be TcG = 0.105 ± 0.003 with the associated 
9cG = 1.2ibl.4 (x^/DOF ~ 1.14 with Q ~ 0.32). The smallness of the error bar associated with 
TcG comes from the fact that Tcross{L) of ^gg/L exhibits a nearly L-independent behavior. 
Large error bar associated with 9cg is merely a consequence of this near L-independence of 
Tcross{L)- Hence, the estimated CG transition temperature, Tcg = 0.105 it 0.003, turns out to 
be higher than the SG transition temperature, Tsg = 0.086 it 0.003, by about 20%, suggesting 
that the spin-chirality decoupling certainly occurs at cr = 0.90. 

One may wonder if the CG susceptibility ratio might be utilized in some way in estimating 
Tcg- As mentioned above, no exact knowledge is available for the value of rjcG- Yet, a finite- 
size scaling analysis performed later in §5 gives an estimate of rjcG — 1-9 at a = 0.9. Thus, we 
show in Fig.5 the temperature dependence of the CG susceptibility ratio xcg/ with 
"HGG = 1-9. The data exhibit a clear crossing at an almost L-independent temperature. A 
power-law extrapolation to L = oo yields an estimate of Tcg = 0.103 it 0.003: See the inset. 
This estimate agrees with the estimate above obtained from ^cg/L, Tcg = 0.105 ± 0.003. 

As can be seen from Fig.3(b), the spin Binder ratio gsG for larger lattices rapidly de- 
creases with increasing L, though a near-merging behavior is observed for smaller lattices at a 
temperature T c± 0.1 close to the CG transition temperature. It should be noticed that, at the 
SG transition temperature, the spin Binder ratio qsg for larger lattices does not exhibit any 
crossing in sharp contrast to the standard crossing behavior observed in the MF SK model or 
the 3D Ising SG, but instead, exhibits only a weak wavy structure somewhat similar to the 
one observed in the 3D SR Heisenberg SG.^^ Presumably, this weak structure would further 
develop into a nontrivial behavior at T = Tsg for still larger lattices. Unfortunately, we cannot 
tell its detailed form at the present stage. 
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In Figs. 6(a) and (b), wc show the chiral-ovcrlap distribution function (a), and the diagonal- 
spin-overlap distribution function (b) at a temperature T = 0.055 which lies well below Tqq 
and TsG- The chiral P{qx) exhibits double peaks at = ±q^^, which tend to diverge with 
increasing L. Unlike the behavior observed in the 3D Heisenberg SG model,^) no central peak 
at = is observed for any size studied. Indeed, with increasing L, the value of P{qx = 0) 
gradually decreases for smaller lattices, but appears to approach a nonzero value of ~ 0.002 for 
larger lattices. Such features of P{qx) are different from the features of the 3D SR Heisenberg 
SG, which exhibits a distinct central peak possibly associated with a one-step-like RSB.^) 
Meanwhile, since the behavior of the Binder ratio of the present model has turned out to be 
entirely different from that of a full RSB,^^^ the ordered state of the present model should not 
simply be regarded as being similar to the one of the MF SK model or of the 3D Ising SG 
model. 

The diagonal-spin-overlap distribution P{qdiag) at a temperature T = 0.055 below Tsg 
shown in Fig.6(b) also exhibits double peaks located at qdiag — ±0.2, which tend to diverge 
with increasing L. The observed diverging peak is the one expected in the SG ordered state 
of the isotropic Heisenberg SG to arise at qdiag = i^q^"^.^^^ Hence, our data of P{qdiag) are 
consistent with a finite SG LR order occurring at this temperature. 

Thus, we have fairly strong numerical evidence of the occurrence of the spin-chirality 
decoupling for a = 0.90, i.e., Tcg = 0.105 ± 0.003 and Tsg = 0.086 ± 0.003. 

4.2 a = 0.85 and a = 0.95 

Next, we investigate the cases of a = 0.85 and 0.95, where the spin-chirality decoupling 
is also likely to occur. The behaviors of the correlation-length ratio and the Binder ratio are 
qualitatively similar to those for a = 0.90 shown above. Hence, we omit exhibiting the cor- 
responding figures just to save space. Meanwhile, a qualitative change occurs in the behavior 
of the SG susceptibility ratio for a = 0.95, where the crossing is no longer observed in the 
temperature and the lattice-size range investigated in contrast to the crossing behavior shown 
in Fig.2. This is demonstrated in Fig.7. 

We show in Fig. 8(a) the size dependence of the crossing temperatures Tcross{L) of ^gg/L, 
^sg/L and XSG / L'^^'^^'-' , for the case of cr = 0.85. The SG transition temperature can be 
estimated from Tcross{L) of Csg/L and of xsG / fairly accurately, each approaching Tsg 
either from above or from below. The combined fit of these two quantities based on eq.(18) 
then yields Tsg = 0.121 ± 0.003 and 9sg = 0.50 ± 0.09 (xVDOF ~ 0.84 with Q ~ 0.52). 
The CG transition temperature is estimated by a power-law fit of TcrossiL) of ^cg/L to be 
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TcG = 0.128 ± 0.003 with Ocg = 0.95 ± 1.58 (xVDOF ~ 0.72 with Q ~ 0.49). Hence, Tsg 
Hes shghtly below Tqg by about 6%. 

In Fig.8(b), we show the size dependence of the crossing temperatures TcrossiL) of £,cg/L 
and of S^sg/L, for the case of cr = 0.95. The SG transition temperature is estimated by a power- 
law fit of TcrossiL) of ^sg/L to be Tsg = OM^^omI, with the associated dsG = 0.36 ± 0.44 
(xVDOF ~ 0.42 with Q ~ 0.52). Note that, for a = 0.95, we have only one type of TcrossiL), 
i.e., that oi£,sG/L, which exhibits a pronounced decreasing tendency with respect to L, leading 
to a rather low estimate of Tsg with the larger error bar. The CG transition temperature is 
estimated by a, power-law fit of Across 

(L) of ^cg/L to be TcG = 0.080 ± 0.005 with the 
associated Ocg = 0.74 ± 0.51 (xV^OF ~ 0.43 with Q ~ 0.65). Hence, Tcg is higher than 
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Fig. 7. The spin-glass susceptibility ratio xsg/L"^ ^^'^ versus the temperature for a = 0.95, with an 
"exact" exponent value r]sG = 3 — 2(t = 1.1. 
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Fig. 8. The (inverse) size dependence of the crossing temperatures of ^cg/^, Csg/-^, and xsg/-^^~''®° 
for a — 0.85 (a) and for a = 0.95 (b). Lines represent power-law fits of the form (18). The CG and 
SG transition temperatures are extrapolated to be Tcg — 0.128 ± 0.003 and TgQ — 0.121 ± 0.003 
for a = 0.85, and Tcg = 0.080 ± 0.005, and Tgc = 0.043to{j4^ for a = 0.95. 



14/28 



J. Phys. Soc. Jpn. 
TsG at a = 0.95. 

4.3 a = 0.8 



Full Paper 




Fig. 9. The correlation-length ratio versus the temperature for the chirality (a), and for the spin (b), 
for (T = 0.8. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. Note 
that the £,cg/L data are given on a semi-logarithmic plot. 




Fig. 10. The spin-glass susceptibility ratio XSG / L"^^^^"^ versus the temperature for a = 0.8, with 
an "exact" exponent value rjsQ = 3 — 2cr = 1.4. The crossing points are indicated by small blue 
arrows. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. 



Next, we study the case of u = 0.8, at which the spin-chirality decoupling ceases to occur 
as we shall see. In Figs. 9-11, we show the correlation-length ratios for the chirality (a) and 
for the spin (b), the SG susceptibility ratio, and the Binder ratios for the chirality (a) and 
for the spin (b), respectively. As can be seen from Fig. 9, at this value of ex, the crossing 
temperatures of the spin S,sg/L come above those of the chiral (,cg/L, the former (the latter) 
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decreases (increases) with increasing L. As can be seen from Fig. 10, the SG susceptibihty 
ratio XsgI^'^~^^'^ shows a clear crossing behavior, where Tcross{L) increases with L. As can 
be seen from Fig. 11, the chiral Binder ratio exhibits a shallow negative dip, while the spin 
Binder ratio exhibits a near- merging behavior below Tsg- 




Fig. 11. The Binder ratio for the chirality (a) and for the spin (b) plotted versus the temperature 
for (7 = 0.8. The red (blue) arrow indicates the bulk chiral-glass (spin-glass) transition point. 
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Fig. 12. The (inverse) size dependence of the crossing temperatures of £,cg/L, ^sg/L and 
XSG I L'^~^^° for (T = 0.8. Lines represent power-law fits of the form (18). The CG and SG tran- 
sition temperatures are extrapolated to be Tcg — 0.158 ± 0.008 and Tsg — 0.159 ± 0.002, i.e., 
TcG = Tsg within the error bar. 



We show in Fig. 12 the size dependence of the crossing temperatures Tcross{L) of ^cg/L, 
^sg/L and xsg/ L'^~^''^^ for the case of cr = 0.8. The SG transition temperature is estimated 
by a combined power-law fit of T^rossiL) of £,sg/L and of XSG / ■ Assuming a common 
Tcross{<y^) = Tsg and 0sg in eq.(18), we get Tsg = 0.158 ± 0.008 with the associated 9sg = 
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0.62 ± 1.95 (x^/DOF ~ 0.43 with Q ~ 0.86). The CG transition temperature is estimated 
via a power-law fit of Tcross{L) of ^cg/L to be Tcg = 0.159 ± 0.002 and Ocg = 0-47 ± 0.11 
(xVDOF ~ 1.14 with Q ~ 0.32). 

44 cr = 0.7 

Next, we study the case of o" = 0.7, which hes close to the lower critical a = |, i.e., the 
boundary between the non-MF and the MF regimes. The behaviors of the correlation-length 
ratios, the SG susceptibility ratio, and the Binder ratios turn out to be more or less similar to 
the ones observed for a = 0.8. One difference is that the spin Binder ratio qsg now exhibits a 
rather clear crossing even for smaller lattices on the positive side of gsG^ as shown in Fig. 13. 




Fig. 13. The Binder ratio for the spin plotted versus the teniperature for a — 0.7. The red (blue) 
arrow indicates the bulk chiral-glass (spin-glass) transition point. 




0.008 



Fig. 14. The (inverse) size dependence of the crossing temperatures of (,cgIL, S.sg/L, and 
XsgI L"^^^^"^ for cr = 0.7. Lines represent power-law fits of the form (18). The CG and SG transition 
temperatures are extrapolated to be Tcg = 0.240 ± 0.002 and Tsg = 0.238 ± 0.002. 
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We show in Fig. 14 the size dependence of the crossing temperatures Tcross{L) of ^sg/L, 
Ccg/I^j and XSG / ■ The SG transition temperature is estimated by a combined power- 
law fit of Tcross{L) of ^sg/L and of XSG / L''^~^^^ ■ Assuming a common Tcrossi^) = Tsg and 
OsG in eq.(18), we get Tsg = 0.238 ± 0.002 and Osg = 0.44 ± 0.21 (xV^OF ~ 0.19 with 
Q ~ 0.94). The CG transition temperature is estimated via a power-law fit of Tcross{L) of 
icG/L to be TcG = 0.240 ± 0.002 with the associated Ocg = 0.28 ± 0.02 (xVDOF ~ 0.0004 
with Q ~ 0.98). Hence, Tsg and Tcg agree within the error bar. 

4.5 cr = 1.0 




Fig. 15. The correlation-length ratio for the chirality (a) and for the spin (b) plotted versus the 
temperature for a ~ 1.0. The red arrow indicates the bulk chiral-glass transition point. Note that 
the S,cgIL data are given on a semi-logarithmic plot. 




Fig. 16. The spin-glass susceptibility ratio XsgI plotted versus the temperature for a = 1.0, 
with an "exact" exponent value rjsG = 3 — 2(t = 1. The red arrow indicates the bulk chiral-glass 
transition point. 
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Next, we turn to the Tsg = regime, i.e., the region of o" > 1. We first study the case of 
a = 1.0, which is just at the upper critical u-value unity, i.e., at the boundary between the 
Tsg > regime and the Tsg = regime. In Figs. 15-17, we show the correlation-length ratios 
for the chirality (a) and for the spin (b), the SG susceptibility ratio and the Binder ratios for 
the chirality (a) and for the spin (b), respectively. As can be seen from Fig. 15, the crossing 
temperatures of both the spin Csg/T and of the chiral Ccg/T curves tend to decrease toward 
lower temperature as L increases. The SG susceptibility ratio xsg/ L'^~^^'^ is shown in Fig. 16. 
Similarly to the one observed for a = 0.95, it does not show a crossing in the temperature and 
the lattice-size range investigated. As can be seen from Fig. 17, the chiral Binder ratio exhibits 
a shallow negative dip. As in the case of o" = 0.90 and 0.95, the chiral Binder ratio of our two 
largest sizes L = 1024 and 2048 exhibits a weak crossing (or a merging) on the positive side 
of gGG around T ~ 0.05, in addition to the one around T ~ 0.08 on the negative side of gcG- 




Fig. 17. The Binder ratio for the chirality (a) and for the spin (b) plotted versus the temperature 
for a — 1.0. The red arrow indicates the bulk chiral-glass transition point. 



We show in Fig. 18 the size dependence of the crossing temperatures Tcross{L) of E,cg/T 
and of ^sg/T. For the spin, the decreasing tendency of Tcross{L) with L becomes pronounced. 
In fact, a power-law fit becomes unstable here, leading to an indefinitely negative T^g"- value. 
Rather, a logarithmic fit of the form expected for the T = transition at the upper-critical 

T,ross{L) = h{\nL + c)-\ (19) 

yields an acceptable fit with 6 ~ 2.1 (x^/DOF= 3.88 and Q = 0.0044) as shown in Fig.l8. This 
observation supports the T = SG transition theoretically expected. For the chirality, a power- 
law fit of Tcross{L) of iGc/L yields Tcg = QMhtomi and Ocg = 0.34 ± 0.34 (xVDOF=0.15 
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and Q = 0.70). Hence, Tcg is likely to be nonzero at o" = 1.0, although the possibility 
of Tcg = cannot be ruled out. Indeed, a logarithmic fit of the form (19) also yields an 
acceptable fit with 6 ~ 0.72 with the x^-value comparable to that of the power-law fit. Note, 
however, even in this case one has Ocg ~ 0.72 << Osg — 2.1, which means that the spin and 
the chirality are decoupled, i.e., ^cg/Csg oo. 

In Fig. 19(a), we show the temperature dependence of the SG correlation length of various 
sizes on a log-log plot together with that of the CG correlation length, which is compared with 
the ones at cr = 1.1 (b) where one also expects Tsg = 0. The slope of an asymptotic straight 
line should give an estimate of the exponent I'sG associated with the T = SG transition. As 
can be seen from Fig. 19(a), significant finite-size effects appear at cr = 1, which prevents us 
from reaching the asymptotic critical regime. Such an unsaturated behavior is also observed 
for the CG correlation length. These behaviors are somewhat in contrast the case of a = 1.1 
shown in Fig. 19(b), where an asymptotic critical behavior seems to be reached both in (^SG 
and ^cG- In the latter cases, one can estimate the exponent usg (or z^cg) associated with the 
T = transition from an asymptotic slope of the data. 

4.6 a = 1.1 

We show in Fig.20 the size dependence of the crossing temperatures Tcross{L) of £,cg/L 
and of £,sg/L for the case of cr = 1.1. Here we find again that the power-law fit becomes 
unstable both for the SG and the CG, yielding indefinitely negative Tsg- or Tec-value. We 
interpret this as suggesting the T = transition lying close to the upper critical cr for both 




Fig. 18. The (inverse) size dependence of the crossing temperatures of S^cg/L, £,sg/L and 
XSG I L'^~^^° for cr = 1.0. Lines for the CG data (red) represent power-law fits of the form (18), 
while lines for the SG data (blue) represent logarithmic-law fits of the form (19). The chiral-glass 
transition temperature is extrapolated to be Tcg = 0-045tQQ27. 
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cases of the SG and the CG. As shown in Fig. 19, the logarithmic plot as employed for a = 1.0 
turns out to work pretty well also for both the SG and the CG at cr = 1.1, with 9 ~ 2.8 for 
the SG, and with 6 ~ 2.4 for the CG. 

In our previous report of Ref.,^^^ we indicated a nonzero Tcg for a = 1.1 on the basis 
of L < 1024 data. In view of an intrinsic difficulty encountered near the upper-critical a in 
distinguishing a Tqg = transition with a large I'cG from a Tqg > transition, however, the 
question of whether the CG transition persists at o" = 1.1 is not completely clear. 

4-7 Phase diagram 

The results obtained in the previous subsections are summarized in the a — T phase 
diagram of Fig. 21. The spin-chirality decoupling occurs in the range 0.8 < o- < 1.1. By 
contrast, the standard spin-chirality coupling behavior with Tsg = Tcg is realized for a ^ 0.8. 
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Fig. 20. The (inverse) size dependence of the crossing temperatures of £,cg/L and of £,sg/L for 
a = 1.1. Lines represent logarithmic-law fits of the form (19). 
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The para-CG phase boundary might go beyond a = 1, touching the T = axis separately from 
the CG-SG phase boundary, although the possibility of it closing just at o" = 1 simultaneously 
with the CG-SG phase boundary cannot be ruled out as shown by the thin dotted line in 
Fig.21. 




2/3 0.7 0.8 0.9 1 1.1 



Fig. 21. The range parameter a versus the temperature T phase diagram of the ID Heisenberg SG 
with a LR power-law interaction decaying with a distance r oc r^'^ . The red (blue) points are 
the chiral TcG (the spin Tsg) transition temperature. The spin-chirality decoupling occurs in the 
range 0.8 < CT < 1.1, while more standard coupling behavior occurs in the range a 0.8. The 
vertical dotted line represents the upper critical a- value, cr — 1. The lower critical cr- value is ct = |. 
The para-CG phase boundary might go beyond a = 1, touching the T = axis separately from 
the CG-SG phase boundary, although the possibility of it closing just at cr = 1 simultaneously 
with the CG-SG phase boundary cannot be ruled out. 



5. Critical properties 

In this section, we study the critical properties of the CG and the SG transitions on the 
basis of a finite-size scaling analysis of our data of the glass susceptibility and the correlation- 
length ratio. Below, we show detailed analyses for the cases of o" = 0.9 and 0.8, each cor- 
responding to the spin-chirality decoupling and coupling regimes, respectively. The critical 
properties for other u-values are also studied, although for the u-values other than a = 0.9 
and 0.8 we quote only the resulting exponent values. 

5.1 a = 0.9 

This value of a corresponds to the spin-chirality decoupling regime. From our analysis 
in the previous section, the CG and the SG transition temperatures were estimated to be 
TcG = 0.105 ± 0.003 and Tsg = 0.086 ± 0.003, respectively. In our following analysis, we fix 
the TcG- and the T^G-values to these best values. The analysis in §4 suggested the presence 
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of a significant correction-to-scaling term. The leading correction-to-scaling exponent was 
estimated to be Ocg = ^CG + 7^ = 1-2 ± 1.4 and 630 = ^^SG + i;^ = O-^^ ± 0.07. In 
our analysis, we will include the effect of the correction-to-scaling by fixing oj + to these 
best values. We then estimate the two independent critical exponents characterizing the CG 
(SG) transitions, i.e., the correlation-length exponent I'cG {^sg) and the critical-point-decay 
exponent r]cG (vsg)- 

We begin with the analysis of the critical properties of the CG transition. We employ the 
finite-size scaling forms for the correlation-length ratio ^cgI^ and for the CG susceptibility 
XCG given by eqs.(13) and (15), respectively. 

As shown in Fig. 22, a reasonably good scaling is obtained both for the CG correlation- 
length ratio £,cg/L and for the CG susceptibility xcG by setting ucG = 4.0 and -qcG = 1-9. 
The associated error bars are estimated by examining by eyes the quality of the fit with 
varying the fitting parameters. We then get vcg = 4.0 it 0.3 and rjcG = 1-9 ± 0.1. 

We perform a similar finite-size scaling analysis also for the spin. We show in Fig. 23 the 
finite-size-scaling plots of the SG correlation-length ratio ^sg/L and of the SG susceptibility 
XSG with the correction term. Both quantities can be scaled reasonably well by setting vsg = 
3.3 and rjsG = 1-2. Our final estimates are then vsg = 3.3 it 0.3 and r]sG = 1-2 ± 0.1. Note 
that the obtained value of rjsG is fully consistent with the analytical expression obtained for 
the LR case rjsG = 3 — 2a. 
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Fig. 22. (Color online) Finite-size-scaling plots of the CG correlation-length ratio £,cg/L (a), and of 
the CG susceptibility xcG (b), for the case of a = 0.9, where the correction-to-scaling effect is 
taken into account. The CG transition temperature and the leading correction-to-scaling exponents 
are fixed to Tcg = 0.105 and ujcg + = 1-2 as determined in §4. The best data collapse for 
^Cg/L is obtained with vcG = 4.0, while that for xcG is obtained with vcg — 4-0 and rjcG = 1-9- 
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Fig. 23. (Color online) Finite-size-scaling plots of the SG correlation-length ratio ^sg/L (a), and of 
the SG susceptibility XSG (b), for the case of cr = 0.9, where the correction-to-scaling effect is taken 
into account. The SG transition temperature and the leading correction-to-scaling exponents are 
fixed to TsG — 0.086 and ujsg + — 0-44 as determined in §4. The best data collapse for ^sg/L 
is obtained with vsg — 3.3, while that for xsG is obtained with vsg — 3.3 and risG — 1-2- 

5.2 a = 0.8 

This value of a corresponds to the spin-chirality coupling regime. From our analysis in 
the previous section, the CG and the SG transition temperatures were estimated to be Tcg = 
0.158 lb 0.008 and Tsg = 0.159 it 0.002, while the leading correction-to-scaling exponent to be 
ecG = ^CG + ^ = 0.62 ± 1.95 and Osg = ^sg + ^ = 0.47 ± 0.11. 

Via the finite-size-scaling analysis shown in Figs. 24 and 25, the CG exponents are deter- 
mined to be vcG = 4.0 lb 0.5 and tigg = 2.0 it 0.1, while the SG exponents are determined 
to be usG = 3.7 ± 0.3 and rjsG = 1-4 ± 0.1. Thus, we get vgg = ^SG within the error bar, 
which is consistent with the expected spin-chirality coupling behavior, i.e., only one diverging 
length scale at the transition. Again, the obtained value of r]sG is fully consistent with the 
analytically obtained expression t]sg = 3 — 2a = 1.4. Note that rjGG and rjsG need not be 
equal even in the spin-chirality coupling case, the spin and the chirality carrying their own 
anomalous dimensions but with a common diverging length scale. 

5.3 Other values of a 

We have performed similar finite-size-scaling analyses for other values of a, and the result- 
ing exponents i^gGi Vcg, ^sg and rfsG are summarized in Fig.26. The analytically obtained 
r/s'G'-value is also included in the figure. As demonstrated in Fig. 19(a), finite-size effects are so 
severe at cr = 1 that we cannot give a reliable estimate of vsg at cr = 1. Nevertheless, as can be 
deduced from the logarithmic fit made in Fig. 18 and from the non-convergent size dependence 
observed in Fig.l9(a), the z/5G-value at cr = 1 could be quite large in the thermodynamic limit. 
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Fig. 24. (Color online) Finite-size-scaling plots of the CG correlation-length ratio ^cg/L (a), and of 
the CG susceptibility XCG (b), for the case of a = 0.8, where the correction-to-scaling effect is 
taken into account. The CG transition temperature and the leading correction-to-scaling exponents 
are fixed to Tea — 0.158 and w + i = 0.62 as determined in §4. The best data collapse for £,cgIL 
is obtained with vcg — 4.0, while that for xcG is obtained with vcg — 4.0 and rjcG — 2.0. 




0.3 -0.2 -0.1 

(T-T.(^)U^SG 



0.3 



Fig. 25. (Color online) Finite-size-scaling plots of the SG correlation-length ratio ^sg/L (a), and of 
the SG susceptibility XSG (b), for the case of cr = 0.8, where the correction-to-scaling effect is taken 
into account. The SG transition temperature and the leading correction-to-scaling exponents are 
fixed to TsG = 0.159 and ujsg + 7^ = 0.47 as determined in §4. The best data collapse for S,sg/L 
is obtained with i/sg — 3.7, while that for xsG is obtained with i^sg — 3.7 and r/sG — 1-4. 



which is not inconsistent with I'sG = c« generically expected at the upper-critical a. 

Several points are to be noticed here, (i) The estimated ijsg agrees well with the analytical 
expression r]sG = 3 — 2a over the entire a < 1 regime, (ii) The estimated r]cG is greater than 
i]SG over the entire a range studied, (iii) At o" = 2/3, usg is expected to approach the MF 
value i^sG = 3. The z^sc-value obtained here toward a = 2/3 appears to be somewhat greater 
than this value. Presumably, a logarithmic correction expected at the lower-critical a = 2/3 
might make an accurate estimate of the exponent difficult around a = 2/3, and the observed 
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exponent might be an effective exponent. Meanwhile, our present estimate of vsG at o" = 0.7, 
^SG = 4.0 lb 0.5, is close to the ones obtained for the ID Ising SG with a LR power-law 
interaction at o" = 0.69, i.e., v = 3.7 it 0.6^^^ or 3.8 it 0.4,^^^ and the one at cr = 0.75, i.e., 
= 4.5 ± 0.2,1*^) 3.3 ± O.4I6) or 4.0 ± 0.5.^^) 
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Fig. 26. (Color online) The a-dependence of the correlation-length exponents vcg and vsg (a), and 
the anomalous-dimension exponents r/cG and rjsc (b). The analytically obtained ?7sG-values, i.e., 
Vsg — 3 — 2(7 for cr < 1 and rjsG = 1 for cr > 1, are also indicated by lines. 



6. Summary and discussion 

We performed a large-scale equilibrium MC simulation on the ID Heisenberg SG with LR 
power-law interactions, paying attention to the SG and the CG orderings and the possible 
spin-chirality decoupling phenomena of the model. This one-dimensional SG model might 
have an advantage over the 3D model that larger linear sizes can be studied. Furthermore, 
by continuously varying and even fine-tuning the power-law exponent a, which plays a role of 
effective "dimensionality", different types of ordering behaviors are realized. 

By calculating various physical quantities including the correlation-length ratio, the sus- 
ceptibility ratio, the Binder ratio, and the overlap distribution function up to the sizes as 
large as L = 4096 for various fi- values in a range of 0.7 < cr < 1.1, we obtained a strong 
numerical evidence for the occurrence of the spin-chirality decoupling behavior in the range 
0.8 S cr S 1.1, while that of the standard spin-chirality coupling behavior with Tsg = Tqg in 
the range cr "S 0.8. Our results are summarized in the temperature-cr phase diagram of Fig. 21. 

Even in the spin-chirality decoupling regime 0.8 < 0- < 1.1, the spin and the chirality often 
behave in a similar way on shorter length scale of L 5 500 ~ 1000, while on longer length 
scale the chirality shows a much stronger ordering tendency than the spin. The observation 
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supports the view of the spin-chirahty couphng behavior at shorter length scale crossing over 
to the spin-chirality decoupling behavior at longer length scale^' 

The Binder ratio exhibits a weak one-step-like RSB feature in the CG ordered state, at 
least in the cr-range where the model exhibits the spin-chirality decoupling behavior. It should 
be noticed, however, that the one-step-like RSB feature observed in the present ID LR model 
is much weaker than the one observed in the CG ordered state of the 3D Heisenberg SG model: 
The dip of the chiral Binder ratio is very shallow and the central peak of the chiral overlap 
distribution function is hardly discernible. By contrast, for a 5 0.8 where the model exhibits 
the standard spin-chirality coupling behavior, we observed a full RSB feature consistently 
with the MF picture. 

We also studied the critical properties of the model based on a finite-size-scaling analysis. 
The resulting exponent values are summarized in Fig. 26. The behavior of the SG exponent 
rjsG is consistent with the analytical result obtained from the RG analysis. 

We try to further examine the possible d-a correspondence. The behavior of ^cg/L of 
the 3D short-range model looks similar to those of cr = 0.9 or cr = 0.95 of the ID LR model. 
To our knowledgeable, no data of the CG and SG correlation length ratios are available 
for higher dimensional Heisenberg SG, which prevents us from making a further comparison. 
Another point to be noticed is that the increasing/decreasing tendency of the size-dependence 
of Tcross{L) of a^cc/L changes from the decreasing (with L) behavior for cr = 1.0 and 0.95 to 
the nearly constant behavior for a = 0.9 and 0.85, and then to the increasing behavior for 
a = 0.8, and 0.7. In comparison with those of the 3D SR models, it appears that the 3D SR 
Heisenberg SG corresponds to a cr-value somewhere between cr = 0.9 or cr = 0.95. Meanwhile, 
the 2D SR Heisenberg SG corresponds to a cr-value around cr = 1.1. 

Overall, the ordering behavior of the Heisenberg SG can roughly be classified into three 
regimes, i.e., the spin-chirality coupling behavior for smaller a (larger d), the spin-chirality 
decoupling behavior for intermediate a (intermediate d) , and the zero-temperature transition 
behavior for larger cr (smaller d). Thus, our present study on the ID LR Heisenberg SG serves 
to provide an overall picture of the ordering behavior of the Heisenberg SG from a wider 
perspective. 
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